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إزالة الصدى والضجيج Echo Cancellation‏ : 
توصيف المشكلة : 
مشكلةل |: (مشكلة تكرار الكلام ف الاتصالات ) 


يف مكالمات الهاتف الدولية والبعيدة المدى بشكل خاص نجد ظاهرة الضجيج ففى حين نجد انه عندما يتكلم الشخص] مع الشخص أ 
نجد ان الشخصل أ يسمع تكرار لصوته وتعزى هذه الظاهرة بسبب تصادم كلا من الاصوات الواردة والصادرة .يستخدم جهاز يدعى 
٩ار"‏ لحل هذه المشكلة والذي يهدف الى فصل اشارات الصوت الصادر عن الصوت الوارد ولكن للأسف هذا الجهاز غير قادر دائما 
على حل هذه المشكلة . 


والشكل التالى يعبر عن نظام لإزالة الضجيج ي دارات الهواتف : 


Voice < E =S +N y__ 
signal; [Delay To 
S | earphone 
noise, n' 
| Adaptive 
circuit filter filter 
To 8 Delay 2 6 
| : ت‎ signal 
earphone 


مشكلة ]: (مشكلة التعرف على صوت الطيار من صوت المحرك ق مركز المراقبة ) 
لنفترض ان لدينا نظام لإزالة الضجيج من الصوت الاصلي التالي : 


Pilot's Volce ۰ 
Pilot's Contaminated with Restored Signal 
Voice Engine Noise 


FN - 


"Error" 


Filtered Noise to Cancel 


Gontaminating Contarnination 


Noise 


Noise Path 


Filter 


Adaptive 
Engine Noise Filter 


Adaptive Filter Adjusts to Minimize Error. 
This removes the engine noiSê from contaminated 
signal, leaving the pilot's voice as the “error.” 
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النظام السابق يفوم بإزالة صوت ضجيج المحرك من الصوت المضفى على صوت فبطان الطائرة .نحن نملك اشارة مفيدة وهي اشارة 
صوت القبطان ي الطائرة ولتكن ۷ .هذه الاشارة تكون موزعة مابين صوت ضجيج المحرك ١‏ »وصوت القبطان وهو يتكلم عبر 
المايكروهون .بكل الاحوال الضجة التى تلوث اشارة صوت الطيار ليست هى نفسها صوت ضجيج المحرك لان الضجة الصادرة عن محرك 
الطائرة ستتخامد بسبب حائط الكبينة وكاتمات اخرى موجودة ف الطائرة . 

امشكلة تكمن ق التخلص من الضجيج لإحتواء صوت الطيار بشكل دفيق . يمكن تحقيق ذلك بتدريب شبكات عصبونية .كي نستطيع 
التنباً باشارة الطيار /محرك ۳ من صوت المحرك ١‏ . ان اشارة صوت المحرك لاتخبر الشبكة العصبية حول صوت الطيار المحتوى ف ٣"‏ 
ولكن الشبكة العصبية يمكن ان تتبناً بصوت المحرك المضاف الى ١‏ . ستحاول الشبكة ان يكون خرجها هو ١‏ ولكن الشبكة تعلم فقط 
صوت الضجة الاصلى (تذكر صوت الضجيج الحقيقى هو صوت المحرك متداخل مع تكتيم حائط الكبينة ١)‏ هذا يعني ان الشبكة 
ستنتج جزء من ١‏ الذي هو مصحح بشكل خطي .هكذا سيكون خطا الشبكة € مساوياً للصوت المشوش (الملوث بالضجيج ) مطروحاً 
منه صوت ضجيج المحرك .ستتعلم الشبكة العصبية حذف صوت ضجيج المحرك بواسطة فلتز مكيف .وفد اثبتت التجارب ان حذف 
الضجيج بواسطة الشبكات العصبية هو افضل من الفلاتر الكلاسيكية لان الضجيج بدل من ان تتم فلترته سيتم الغائه . 


يمكن تفسيم المشكله التى نحنا بصدد حلها الى عدة مراحل : 


تعريف اشارة الصوت اج١‏ عاS‏ عVoic‏ 


تعريف اشارة الضجيج ١2|‏ عاك عءأهل (ضجيج ابتدائي ) 


| تعريف فلتر مستخدم لإنشاء ضجة ملوثة مءزمل ١ع‏ اه٤‏ (ضجة حقيقة ) 


Contaminate Signal ڈوlم تعریف إشارة‎ 


تعريف شبكة عصبونية ولتكن ع ال۸ (كونها الاشهر في مجال معالجة الإشارة ) 


إظهار النتائج بشكل بياني 
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المرحلةا | ١‏ تعريف اشارة الصوت اة عا؟S‏ عVoic‏ 


اول دخل يجب تعريفه هو صوت الطيار / إن اشارة صوت الطيار يمكن تحديدها وفق المطال والتردد Frequency & a p|i†u d€‏ إj‏ 
تحدید کلاً من المطال والتردد لصوت الإنسان هو هام جدا فالمطال يحد مدى علو الصوت وانخفاضه 0655 |0٥ U‏ بينما يحدد التردد درجة الصوت 


) اأموتعطى معادلة صوت الطيار (الإنسان ) وفق ما يلي (المعادلة وفق منهج ومقرر الاتصالات الرقمية للسنة الثانية وضعها العالم بابلينسكي‎ ٣ 


V(t) = (1 + 0.2 sin(Vamp) t). sintfMS(1 + 0.2 cos(Vfreq t))t) 4 
. دلالات المعادلة‎ 


الزاوی حامل الاشارة ۲ ۲۲1€ Angular Freq of the C3‏ .بحیث تکون قیمته مساویة ل p۴‏ 2 


ان التردد الزاوي لحامل الإشارة هو التردد الاساسي لاشارة الكلام أ2ة ١ع‏ أك 50886٤1‏ . إن تردد الإشارة ؟ تم اختيارها لتكون 
KHzİ Î‏ 


Vfreq =2 pm 


دل ت دد ال mn‏ ّ عه م" ل¿ تردد الإأشارة . مه E.‏ 
هي معادلة تردد التعديل » م هي عملياً مجموعة من القيم الاعتباطية اقل من تردد الإشارة . وتعطى قيمة 300112 ف هذ 
اللالة من دراستنا. 


Vamp=2Pfa 
. 200۸7 هي معادلة تردد مطال التعديل »ان مطال التعديل دأ ياخذ قيم عديدة اقل من تردد الاشارة الاصلي . نعطيه قيمة‎ 
=٣ ان معادلة الزمن هي‎ 
. ه5 ويكون آ [] ميلي ثانية‎ م١8‎ [١€ بحيث ان ك هي زمن التقطيع‎ 
: باستعمال هذه المتحولات يمكن تعريف اشارة صوت ف ماتلاب والكود التالي لتطبيق هذا الكلام‎ 


f=4e3; % set a signal freq of 4 khz 
fm=300; %set freq modulation of 300 hz 
fa=200; %set an amplitude modulation of 200 hz 


ts=2e-5 ; %set a sampling time of 0.2 mele sec 

N=400; % set 400 sampling points 

t=(O:N-1)*ts ; %set a discrete time from O to 10 mele sec 
v=(1+.2*sin(2*pi*fa*t)).*sin(2*pi*f*(1+.2cos(2*pi*fm*t)).*t); voice signal 
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المرحلةا | ٭تعريف اشارة ضıج Noise Signal‏ : 


ان اشارة الضجيج ۸ الآن يمكن لها ان تعرف بأنها اشارة متلثية للتردد 1ikhz= fh‏ المدمج مع عناصر من ضوضاء عشوائية مثلا 
(()engthا,1) ran dn‏ ان اشارات اسنان المنشار كما يطلق عليها هي اشارة تكرارية وتم اختيارها لسهولة الدراسة .ان نوع الإشارات لس ها 


شكل اشارة اسنان المنشار 00 5W‏ قبل وبعد الضجيج .(المصدر موقع (math WOKS‏ 


إن اشارة اسنان المنشار سيكون ها قمم عند الأ ] السالبة والموجبة ‏ وق تكرار دوري لكل 2 .سيتم إضفاء ضجيج عشواني هذه 
الإشارة كما ف الشكل الموضح آنفاً . هذا الامر سيزودنا باشارة ضجيج متغيرة باستمرار تتقلب للأعلى والادنى .هذه الإشارة هي 
مناسبة جداً لاجراء اختبار ١€‏ |303 بسبب ان الطيار سيتكلم بصوت اخفض من صوت هدير المحرك . 


إن هذا الخليط بين تكرار الإشارة والإشارة العشوائية ت اختياره بسبب ان الصوت المراد معالحته هو صوت ناتج لشيئين اثنين : 
صوت الطيار (صوت بشري ) وصوت هدير المحرك الذي يعمل بشكل تكراري ودوري لذلك كانت تلك التكرارية قى اشارة اسنان 
المنشار ماهي إلا لتسهيل دراسة ضجيج المحرك . 


المرحلةا | و تعريف فلتر مستخدم لإنشاء ضجة ملوثة موز dم۲هام٣‏ (ضجة حقيقة ) 


إن الضجيج المكتوم الذي يسمعه الطيار يمكن تكوينه بتلوين الضجيج باستخدام فلتر خطي ۴|۴ الذي يملك شعاع مجهول ج في الحالات الحفيقة 
يكون الشعاع ج مجهولاً ولكن في حالة دراستنا فإنه سيكون معروف القيم وذلك بسبب اختيار فلتر ٣ا۴‏ » هذا الامر يعني ان اشارة الضجيج الضعيفة 


تحتاج الى تنبوءات من قبل شبكة "٥‏ اة هه . ان فلتر ٣ا۴‏ يستخدم بسبب ان الضجيج يسافر من خلال بيئة غير خطية قبل ان تصل نسخة خافتة الى 
ادن الطيار .ان فلتر ٣ا۴‏ هو ابسط نماد ج الفلاتر من اجل طریق انتقال الضجيج 1 


fn=1e3 ; %set signal frequency of 1Khz 
n=randn(1,length(t))+sawtooth(2*pi*fn*t,0.7); %create noise 


a=[1 -0.6 -0.3|]; %filter coefficients 
c=filter(a,1,n); create output signal c from Fir filter 
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المرحلةا | تعریف إشارة ذوûڎة Contaminate Signal‏ 
الكود التالي يبين كيف يمكن للإشارة الملوثة ع ان يتم اضافتها الى صوت الطيار لننشىء اشارة ملوثة حقيقة 


m=v+c; %pilots voice contaminated with engine noise 


يمكن رسم كلا من الإشارتين صوت الطيار وصوت ضجيج المحرك وفق تعليمة اهام لكل من ۷,م وبالتالي ينتج لدينا المخطط التالي : 


o CE engine sg nal‏ —ے 
voice signal‏ _—— 


0 1 2 3 4 5 6 7 8 


timêelmês) ¥ 1 


: :كيفية تعامل ۹|۹6 مع ملفات الصوت‎ ]mPاemen‎ in 
: ”تحمیل ملف صوت الى برنامج طaاةص :يقوم البرنامج بتحميل ملفات الصوت وتحويلها الى مصفوفة ×أ٣ ةج" وفق التعليمة التالية‎ 
[sound,fs]=wavread('c:\test.wav'); %read a wav file from a file 


في حال لم نقم بوضع ; "فاصلة منقوطة " سنجد اننا قمنا بتوليد مصفوفة ارقام تمثل ترددات الموجات الصوتية للملف »ووضع تلك القيم في 
المصفوفة المنشئة حديثاً ل" هء في حال قمنا باستعراض المصفوفة ل مء سنجد الارقام المكونة لملف الصوت . ويمكن تخزين جميع متحولات 
البرنامج المستخدمة في ملف ذو لاحفَة ج ص. وفق التعليمة التالية لمشاهدة قيم المتحول 4ء في ملف 


filePath ='c:\TestFilelnfo.mat'; %the path of the wavfile plus the file name 
save(filePath,'sound'); %save the sound wav file into mat file 
5 ک“ ان نشاهد حجم | 1 وفق اأ أب‎ ٤ 
size(varName) 

Left=sound(:,1); % sound Matrix contain 2 columns for left and right channel 

لنقم برسم ترددات الصوت بدلالة الزمن وفق التعليمات التالية (الصورة الثانية للإيضاح اكثر ) 
سڪ time=(1/44100)*length (left); a‏ 

USE HS | kA 0 DBD O 4 - 3| O E | a 


1 Note new toolbar buttons: data brushing & linked plats a al Play widea 


t=linspace(0,time,length(left)); 
plot(t,left) 


xlabel('time (sec)'); 


relative signal strength 


ylabel('relative signal strength’) 
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time=(1/44100)*2000; J Figure 1 


File Edit Wiew Insert Tools Desktop Window Help 
USE Rة‎ | kk | ك‎ & 3 DBD O £ || O E | a 
O Note new toolbar buttons: data brushing & linked plots a al Play video 


t=linspace(0,time,2000); 0.2 
plot(t,left(1:2000)) 


xlabel('time (sec)'); 


relative signal strength 


ylabel('relative signal strength’) 


0-005 001 0015 002 0025 003 0035 004 0045 005 
time (sec) 


يمكن سماع الصوت وفق قناة صوت واحدة ولتكن السار دة يكون ذلك غن طريق تعامة 
soundsc(left,fs)‏ 


يمكن اضافة العديد من انوع الضجيج وتشوهات الصوت والضجيج التالي يبين تكرار للصوت كما في مكالمات الهاتف الدولية : 


eftout=left; % set up a new array, same size as old one 


N=10000; % delay amount N/44100 seconds 
for n=N+1:length (left) 


leftout(n)=left(n)+left(n-N); % about % second echo 


ویمکن المقارنة بين ملفي الصوت الاصلي والمشوه : 


>> soundsc(left,fs) 
>> soundsc(leftout,fs) 


المرحلةا | تعريف شبكة عصبونية ع iاةكA:‏ 


بسبب ان بارميترات ضجيج الفلتر ناهيك عن كون ج مجهولة فإن من الضروري استنتاجه لتحقيق عملية إلغاء الضجيج .ان البراميتر يمكن 
استنتاجه بواسطة شبكة ,٥‏ اج له .فيما بعد فإن الضجيج المسبب لتشويه صوت الطيار يمكن استنتاجه . في هذه الحالة فإن شبكة "٥‏ iاة‏ له يمكن 
تشكيلها بواسطة تعريف شبكة ذات خرج وحيد ولتكن عم . بينما قيم الدخل هي مجموعة تتراوح بين اكبر قيمة لضجيج الإشارة الى اصغرها . هذا 
الامر بسبب ان اشارة الضجيج تدخل الى شبكة ع iاةلa‏ . سيتم وضع 3 وحدات تاخير مع اوزانها الخاصة التي تعرف من اجJ tapped delay‏ 
مما . من اجل اشارة اسنان المنشار (الضجيج ) فإنه من الكافي ان نملك نيرون ۲٥٣‏ ں مل مع 4 اوزان . تم إيجاد هذا الامر وفق التجريب 
ومشاهدة الاخطاء خلال التجريب لذلك فإن خط التاخير سيتصل مع مصفوفة اوزان الشبكة من خلال تاخير 0و 1 و 2و 3 وحدات زمن .ان 
مصفوفة الخلايا ام تعرف قيم الخرج الابتدائي للتاخير ومن المهم التاكد على ان معامل الانزياح ءهأط يملك قيم هي عبارة عن مجموعة من القيم 
المرغوبة خصوصاً صفر . 
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المر : ل | تنفبد لية تدریب ال رک 0 


ان دخل الشبكة هو م والخرج المطلوب هو م يحتاج الى ان نعرفه كعناصر من مصفوفة خلايا رج٣٣ج‏ اام . ان مصفوفة الخلايا هي بنى 
E E‏ ق کر اک ع ا ا ا ا 
هدفه عع ۲ه من العمود ¡ في خلايا الهدف . 


من اجل دخل معطى من ع٣ا‏ رمامك ١عممها‏ اي 3-أ1,"1-2,۸-أ,ام الشبكة ستتدرب لتعطي خر ج ا" 


بتحديد 1 كقيمة لسلسة الدخل فان الشبكة المعرفة يمكن تدرببها بد٤‏ من شرط التاخير المعرف آنفاً . 


net=newlin([min(n),max(n)],1); %assume input values 
net.inputWeights{1,1}.delays=[0 1 2 3 |]; %3 delays 
net.IW{1,1}=[1 1 1 1|; %give the various weights these vals 
net.b{1}=[0]; %set bias to O 


Pi={1 1 1} %define the initial output of the delays 
n=num2cell(n) %change arrays into the required cells 
m=num2cell(m); 

net.adaptparam.passes=1; %specify 1 pass 

[net,y,E pf,af]=adapt(net,n,m,pi); do training 


يمكن حساب الاخطاء من التالي : 


temp=zeros(400,1); 
for i=1 :400 
temp(i,1)=y{1,i}(1,1); 
end 
temp=temp’; %transpose the array 
m=cat(2,m{:}); 
error= (v-(m-temp)); %calculate error of signal 


error2=0; 

temp=zeros(400,1); 

for i=1:400 
error2=error2+abs(error(i)) 

end 

error3=error2/400 


إذا صغرنا الثابت بشكل كبير لن تتغير الإشارة الناتجة عن الطرح متل 0.0001 


اذا كبرتا متل 0.6 و أكتر بيصير في تشويش فافضل قيمة هي ما بين 0.01 و ال 0.35 تقريبا 
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3- بعد إجراء التدريب نجد ان الضجيج نم إلغاءه وفق المخططات التالية : 


المخطط التالي يبين معدل الاخطاء : 


[“ ااا ا ا aE‏ 


0 0 1 ل 1 ا‎ N 0 0 N 1 


YT 


ime(ms) و‎ 


المخطط التالى يبين الضجيح الملوث وخرج شبكة 2|1۸€ 3d‏ 


LL | 
—— coloured noisê 
—_ oDltput of Adalinê 


ime(ms) 3 


المخطط التالى ييبين حالة اشارة الصوت واشارة الصوت الملوثة . 


ime(ms) 


9ع P2‏ | محمد وسيم أبوزينة 


